Generation of Hyperentangled Photon Pairs 
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We experimentally demonstrate the first quantum system entangled in every degree of freedom 
(hyperentangled). Using pairs of photons produced in spontaneous parametric down- conversion, we 
verify entanglement by observing a Bell-type inequality violation in each degree of freedom: polar- 
ization, spatial mode and time-energy. We also produce and characterize maximally hyperentangled 
states and novel states simultaneously exhibiting both quantum and classical correlations. Finally, 
we report the tomography of a 2x2x3x3 system (36-dimensional Hilbert space), which we believe is 
the first reported photonic entangled system of this size to be so characterized. 
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Entanglement, the quintessential quantum mechanical 
correlations that can exist between quantum systems, 
plays a critical role in many important applications in 
quantum information processing, including the revolu- 
tionary one-way quantum computer quantum cryp- 
tography Q, dense coding |>] and teleportation As 
a result, the ability to create, control and manipulate 
entanglement has been a defining experimental goal in 
recent years. Higher-order entanglement has been real- 
ized in multi-particle Q and multi-dimensional 0,0,11,0 
systems. Furthermore, two-component quantum systems 
can be entangled in every degree of freedom (DOF), or 
hyperentangled [l0[. These systems enable the imple- 
mentation of 100%-efflcient complete Bell-state analysis 
with only linear elements [TTf] and techniques for state 
purification [l2[. In addition, hyperentanglement can 
also be interpreted as entanglement between two higher- 
dimensional quantum systems, offering significant advan- 
tages in quantum communication protocolsfe.g., secure 
superdense coding [13] and cryptography [~L4l]). 

Photon pairs produced via the nonlinear optical pro- 
cess of spontaneous parametric down-conversion have 
many accessible DOF which can be exploited for the pro- 
duction of entanglement. This was first demonstrated 
using polarization [Til, [l6[ , but the list expanded rapidly 
to include momentum (linear [17], orbital [y] , and trans- 
vers e Il8| spatial modes), energy-time [19j and time- 
bin [20[ , simultaneous polarization and energy-time [2l( , 
and recently, simultaneous polarization and 2-level linear 
momentum [22|. In this work, we produce and charac- 
terize pairs of single photons simultaneously entangled in 
every DOF -polarization, spatial mode and energy-time. 
As observed previously Q, photon pairs from a single 
nonlinear crystal are entangled in orbital angular momen- 
tum (OAM). Moreover, polarization entangled states can 
be created by coherently pumping two contiguous thin 
crystals [23|, provided the spatial modes emitted from 
each crystal are indistinguishable. Finally, the pump dis- 
tributes energy to the daughter photons in many ways, 
entangling each pair in energy; equivalent ly, each pair 



is coherently emitted over a range of times (within the 
coherence of the continuous wave pump). We show our 
two-crystal source can generate a2x2x3x3x2x2- 
dimensional hyperentangled state [10[, approximately 

(\HH) + \VV)) ® (\rl) + a\gg) + \lr)) ® (\ss) + \ff)) . 
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Here H (V) represents the horizontal (vertical) photon 
polarization; |Z), \g) and |r), represent the paraxial spa- 
tial modes (Laguerre-Gauss) carrying —ft, 0, and +ft 
OAM, respectively [24[; a describes the OAM spatial 
mode balance prescribed by the source [25] and selected 
via the mode-matching conditions; and \s) and |/), re- 
spectively, represent the relative early and late emission 
times of a pair of energy anticorrelated photons [191 ]. 

The most common maximally entangled states are 
the 2-qubit Bell states: $± = (|00) ± |ll))A/2 and 
^ ± = (|01) ± |10))/V2, in the logical basis |0) and |1). 
By collecting only the ±ft OAM state of the spatial sub- 
space, the state (pp) becomes a tensor product of three 
Bell states $p oln (8> (8> 3>tt e - As a preliminary test of 
the hyperentanglement, we characterized the polarization 
and spatial mode subspaces by measuring the entangle- 
ment (characterized by tangle T 0), the mixture (char- 
acterized by linear entropy Sl(p) = f [1 - Tr(p 2 )] [27|). 
and the fidelity F(p, pt) = (Tr(^/ \/ptp^/pt)) 2 of the mea- 
sured state p with the target state p t = \ipt)(ipt\- We 
consistently measured high-quality states with tangles, 
linear entropies, and fidelities with <I> + of T = 0.99(1), 
S l = 0.01(1) and F = 0.99(1) for polarization; and 
T = 0.96(1), S l = 0.03(1) and F = 0.95(1) for spatial 
mode, significantly higher than earlier results [18]. The 
experiment is illustrated in Fig. 1. A 120-mW 351-nm 
Ar + laser pumps two contiguous /3-barium borate (BBO) 
nonlinear crystals with optic axes aligned in perpendic- 
ular planes [23[. Each 0.6-mm thick crystal is phase- 
matched to produce Type-I degenerate photons at 702 
nm into a cone of 3.0° half-opening angle. The first (sec- 
ond) crystal produce pairs of horizontally (vertically) po- 



2 



larized photons and these two possible down- conversion 
processes are coherent, provided the spatial modes emit- 
ted from each crystal are indistinguishable. With the 
pump focused to a waist at the crystals, this constraint 
can be satisfied by using thin crystals and "large" beam 
waists (large relative to the mismatch in the overlap of 
the down-conversion cones from each crystal [23[). How- 
ever, the OAM-entanglement is maximized by balancing 
the relative populations of the low-valued OAM eigen- 
states [25j, which requires smaller beam waists to im- 
age a large area of the down- conversion cones. Here 
we compromise by employing an intermediate waist size 
(~ 90 fim) at the crystal. Mode-matching lenses are 
then used to optimize the coupling of the rapidly diverg- 
ing down-conversion modes into single-mode collection 
fibers. 

The measurement process consists of three stages of 
local state projection, one for each DOF. At each stage, 
the target state is transformed into a state that can 
be discriminated from the other states with high ac- 
curacy. Specifically, computer-generated phase holo- 
grams [28| transform the target spatial mode into the 
pure gaussian (or 0-OAM) mode, which is then filtered 
by the single-mode fiber @ (Fig. lb). After a po- 
larization controller to compensate for the fiber bire- 
fringence, wave plates transform the target polarization 
state into horizontal, which is filtered by a polarizer 
(Fig. Id). The analysis of the energy-time state is re- 
alized by a Franson-type [19] polarization interferome- 
ter without detection post-selection [2l[. The matched 
unbalanced interferometers give each photon a fast |/) 
and slow \s) route to its detector. Our interferometers 
consisted of L ~ 11-mm quartz birefringent elements, 
which longitudinally separated the horizontal and ver- 
tical polarization components by An quartz L ~ 100 /im, 
more than the single-photon coherence length (A 2 /AA^ 
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FIG. 1: Experimental setup for the creation and analysis of 
hyperentangled photons, (a) The photons, produced using 
adjacent nonlinear crystals (BBO), pass through a state fil- 
tration process for each DOF before coincidence detection. 
The measurement insets show the filtration processes as a 
transformation of the target state (dashed box) and a filter- 
ing step to discard the other components of the state (dot- 
ted box), (b) Spatial Filtration (spa): hologram (holo) and 
single- mode fiber (smf). (c) Energy-Time Transformation 
(e-t): thick quartz decoherer (dec) and liquid crystal (LC). 
(d) Polarization Filtration (poln): quarter- wave plate (qwp), 
half- wave plate (hwp) and polarizer (pol). 



50 /am with AA = 10 nm from the interference filters) 
but much less than the pump-photon coherence length 
(~ 10 cm). We rely on the photons' polarization entan- 
glement \HH) + \VV) to thus project onto a two-time 
state (\Hs,Hs) -\- e^ 6l+s ^ \V f , Vf)), where 5 X and S 2 are 
controlled by birefringent elements (liquid crystals and 
quarter- wave plates) in the path of each photon [2l| . Fi- 
nally, by analyzing the polarization in the ±45° basis, 
we erase the distinguishing polarization labels and can 
directly measure the coherence between the \ss) and \ ff) 
terms, arising from the energy-time entanglement. 

To verify quantum mechanical correlations, we 
tested every DOF against a Clauser-Horne-Shimony-Holt 
(CHSH) Bell inequality @. The CHSH inequality places 
constraints (S < 2) on the value of the Bell parameter S, 
a combination of four two-particle correlation probabili- 
ties using two possible analysis settings for each photon. 
If S > 2, no separable quantum system (or local hid- 
den variable theory) can explain the correlations; in this 
sense, a Bell inequality acts as an "entanglement wit- 
ness" [30|. To measure the strongest violation for the 
polarization and spatial-mode DOFs, we determined the 
optimal measurement settings by first tomographically 
reconstructing the 2-qubit subspace of interest; we em- 
ploy a maximum likelihood technique to identify the den- 
sity matrix most consistent with the data [27[. 

Table [J shows the Bell parameters measured for the po- 
larization, spatial mode, and energy-time subspaces, with 
various projections in the complementary DOF. We see 
that for every subspace, the Bell parameter exceeded the 
classical limit of S = 2 by more than 20 standard devia- 
tions (a), verifying the hyperentanglement. For both the 
polarization and spatial-mode measurements, we traced 
over the energy-time DOF by not projecting in this sub- 
space. We measured the polarization correlations while 
projecting the spatial modes into the orthogonal basis 
states (|Z), and |r)), as well as the superpositions 
\h) = (|0 + \r))/y/2 and \v) = - \r))/y/2). The mea- 
sured Bell parameters agreed (within ~ 2 a) with predic- 
tions from tomographic reconstruction and violated the 



TABLE I: Bell parameter S showing CHSH-Bell inequality 
violations in every degree of freedom. The local realistic limit 
(S < 2) is violated by the number of standard deviations 
shown in brackets, determined by counting statistics. 
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inequality by more than 30 a. In the spatial mode DOF, 
the correlations for the state a were close to maximal 
(S = 2^/2 « 2.83), also in agreement with predictions 
from the measured state density matrix. In addition, 
we tested Bell inequalities for non-maximally entangled 
states in the OAM-subspace: a\gg) + \rl) and a\gg) + \lr); 
the measured Bell parameters in this case were slightly 
smaller (5%, max.) than predictions from tomographic 
reconstruction [3l|, yet still 20 a above the classical limit. 
Finally, our measured Bell violation for the energy-time 
DOF using particular phase settings is in good agreement 
with the prediction (S = 2\^V) from the measured 2- 
photon interference visibility V = 0.985(2). 

The polarization and spatial-mode state was fully char- 
acterized via tomography [27]. We performed the 1296 
linearly independent state projections required for a full 
reconstruction in the (2 <g) 3) <g) (2 <g> 3) Hilbert space con- 
sisting of two polarization and three OAM modes for each 
photon. The measured state (Fig. 2) overlaps the antic- 
ipated state (polarization and spatial DOFs of Eq. 1) 
with a fidelity of 0.69(1) for a = 1.88e°- 16i7r (numerically 
fitted), and Sl = 0.46(1), suggesting the difference arises 
mostly from mixture. Treating the photon pairs as a 
six-level two-particle system, we can quantify the entan- 
glement using the negativity N [33[ . In this 6 6 Hilbert 




FIG. 2: (color online). Measured density matrix (p) and close 
pure state (|^ p ) ~ $+ oln ® (\lr) + a\gg) + \rl)) with a = 
1.88e°- 16i7r ) of a (2x2x3x3)-dimensional state of 2-photon 
polarization and spatial mode [321 ]. 
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F(p, p t )= 0.974(1) 

N(p) = 2.910(3) 
N( Pt ) = 3 
S L (p) = 0.039(2) 
S L ( Pt ) =0 
Subspaces: 
T poln = 0.945(2) 
T spa = 0.943(3) 
S L , pain = 0.035(2) 
S L , spa = 0.033(2) 
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F(p, Pt )= 0.955(2) 

N(p) = 2.884(9) 
N(p t ) =3 
S L (p) = 0.048(5) 
S L (pt) =0 
Subspaces: 
Tp i n = 0.925(6) 
T spa = 0.937(4) 
Si,, pom = 0.048(4) 
S L , spa = 0.034(3) 
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p t )= 0.972(1) 

N(p) = 0.983(4) 
N( Pt ) =1 
S L (p) = 0.549(1) 
S L (pt) =0.533 
Subspaces: 
T P om = 0.000(1) 
T spa = 0.938(3) 
S L , pain = 0.675(1) 
S L , spa = 0.035(2) 
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F(p, p t )= 0.979(1) 

N(p) = 0.969(2) 
N(p t ) =1 
S L (p) = 0.804(1) 
S L (pt) =0.800 
Subspaces: 
Tpoln = 0.000(1) 
T spa = 0.930(3) 
Si,, pom = 0.998(1) 
S L , sp* = 0.039(2) 



FIG. 3: (color online). Measured density matrices (real parts) 
of (2x2x2x2)-dimensional states of 2-photon polarization 
and (+1, — l)-qubit OAM [33]. For each state, we list: the 
target state p t , the fidelity F(p,p t ) of the measured state p 
with the target pt, their negativities and linear entropies, and 
the tangle and linear entropy for each subspace. The neg- 
ativity for two-qubit states is the square root of the tangle. 
The magnitudes of all imaginary elements, not shown, are less 
than 0.03. 
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space, N ranges from (for separable states) to 5 (for 
maximally-entangled states), and the fitted state above 
has N & 4.44 Our measured partially mixed state has 
N = 2.96(4), indicating strong entanglement. The spatial 
mode alone has N= 1.14(2), greater than the maximum 
(N = 1) of any two-qubit system. Thus, our large state 
possesses 2-qubit and 2-qutrit entanglement. 

We also selected a state (neglecting the \gg) compo- 
nent, Fig. 3a) maximally entangled in both polarization 
and spatial mode, that had F = 0.974(1) with the target 
^poin®^spa- By tracing over polarization (spatial mode), 
we look at the measured state in the spatial mode (po- 
larization) subspaces. The reduced states in both DOFs 
are pure (Sl< 0.04) and highly entangled (T>0.94). 

With this precise source of hyperentanglement, we 
have the flexibility to prepare nearly arbitrary polar- 
ization states [34], and to select arbitrary spatial- mode 
encodings. For example, we also generated a different 
maximally entangled state: ^p oln (Fig. 3b). By 
coupling to and tracing over the energy-time DOF us- 
ing quartz decoherers [341 ]. we can add mixture to the 
polarization subspace, allowing us to prepare a previ- 
ously unrealized state that simultaneously displays clas- 
sical correlations in polarization and maximal quan- 
tum correlations between spatial modes (Fig. 3c): p « 
\{\HH){HH\ + \VV){VV\) ® |$+ pa )(*+al- We were 
also able to accurately prepare the state p t = jI po \ n ® 
l^spaK^spal? with no polarization correlations at all (i.e., 
completely mixed or unpolarized), while still maintaining 
near maximal entanglement in the spatial DOF (Fig. 3d). 



We report the first realization of hyperentanglement of 
a pair of single photons. The entanglement in each DOF 
is demonstrated by violations of CHSH-Bell inequalities 
of greater than 20a. Also, using tomography we fully 
characterize a 2 (g) 2 <g) 3 (g) 3 state, the largest quantum 
system to date. In restricted (2 x 2 x 2 x 2)-dimensional 
subspace, we prepare a range of target states with un- 
precedented fidelities for quantum systems of this size, 
including novel states with a controllable degree of cor- 
relation in the polarization subspace. These hyperentan- 
gled states enable 100%-efflcient Bell-state analysis [111 ] . 
which is important for a variety of quantum information 
protocols |3, LL3[ . Because the spatial mode and energy- 
time DOFs are infinite in size, we envision examining 
even larger subspaces, encoding higher-dimensional qu- 
dits [7, 8]. Finally, we note that the pair wise mechanism 
of the down-conversion process inherently produces 
entanglement in photon number [35]. 
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